Increasing the power of genome wide association studies in natural populations using repeated measures : evaluation and implementation

1. Genomewide association studies (GWAS) enable detailed dissections of the genetic basis for organisms' ability to adapt to a changing environment. In long-term studies of natural populations, individuals are often marked at one point in their life and then repeatedly recaptured. It is therefore essential that a method for GWAS includes the process of repeated sampling. In a GWAS, the effects of thousands of single-nucleotide polymorphisms (SNPs) need to be fitted and any model development is constrained by the computational requirements. A method is therefore required that can fit a highly hierarchical model and at the same time is computationally fast enough to be useful. 2. Our method fits fixed SNP effects in a linear mixed model that can include both random polygenic effects and permanent environmental effects. In this way, the model can correct for population structure and model repeated measures. The covariance structure of the linear mixed model is first estimated and subsequently used in a generalized least squares setting to fit the SNP effects. The method was evaluated in a simulation study based on observed genotypes from a long-term study of collared flycatchers in Sweden. 3. The method we present here was successful in estimating permanent environmental effects from simulated repeated measures data. Additionally, we found that especially for variable phenotypes having large variation between years, the repeated measurements model has a substantial increase in power compared to a model using average phenotypes as a response. 4. The method is available in the R package RepeatABEL. It increases the power in GWAS having repeated measures, especially for long-term studies of natural populations, and the R implementation is expected to facilitate modelling of longitudinal data for studies of both animal and human populations.

Association of soluble tumor necrosis factor receptors 1 and 2 with nephropathy, cardiovascular events, and total mortality in type 2 diabetes

AIMS/HYPOTHESIS: Soluble tumor necrosis factor receptors 1 and 2 (sTNFR1 and sTNFR2) contribute to experimental diabetic kidney disease, a condition with substantially increased cardiovascular risk when present in patients. Therefore, we aimed to explore the levels of sTNFRs, and their association with prevalent kidney disease, incident cardiovascular disease, and risk of mortality independently of baseline kidney function and microalbuminuria in a cohort of patients with type 2 diabetes. In pre-defined secondary analyses we also investigated whether the sTNFRs predict adverse outcome in the absence of diabetic kidney disease. METHODS: The CARDIPP study, a cohort study of 607 diabetes patients [mean age 61 years, 44 % women, 45 cardiovascular events (fatal/non-fatal myocardial infarction or stroke) and 44 deaths during follow-up (mean 7.6 years)] was used. RESULTS: Higher sTNFR1 and sTNFR2 were associated with higher odds of prevalent kidney disease [odd ratio (OR) per standard deviation (SD) increase 1.60, 95 % confidence interval (CI) 1.32-1.93, p < 0.001 and OR 1.54, 95 % CI 1.21-1.97, p = 0.001, respectively]. In Cox regression models adjusting for age, sex, glomerular filtration rate and urinary albumin/creatinine ratio, higher sTNFR1 and sTNFR2 predicted incident cardiovascular events [hazard ratio (HR) per SD increase, 1.66, 95 % CI 1.29-2.174, p < 0.001 and HR 1.47, 95 % CI 1.13-1.91, p = 0.004, respectively]. Results were similar in separate models with adjustments for inflammatory markers, HbA1c, or established cardiovascular risk factors, or when participants with diabetic kidney disease at baseline were excluded (p < 0.01 for all). Both sTNFRs were associated with mortality. CONCLUSIONS/INTERPRETATIONS: Higher circulating sTNFR1 and sTNFR2 are associated with diabetic kidney disease, and predicts incident cardiovascular disease and mortality independently of microalbuminuria and kidney function, even in those without kidney disease. Our findings support the clinical utility of sTNFRs as prognostic markers in type 2 diabetes.

Genome-wide association mapping in a wild avian population identifies a link between genetic and phenotypic variation in a life-history trait

Understanding the genetic basis of traits involved in adaptation is a major challenge in evolutionary biology but remains poorly understood. Here, we use genome-wide association mapping using a custom 50 k single nucleotide polymorphism (SNP) array in a natural population of collared flycatchers to examine the genetic basis of clutch size, an important life-history trait in many animal species. We found evidence for an association on chromosome 18 where one SNP significant at the genome-wide level explained 3.9% of the phenotypic variance. We also detected two suggestive quantitative trait loci (QTLs) on chromosomes 9 and 26. Fitness differences among genotypes were generally weak and not significant, although there was some indication of a sex-by-genotype interaction for lifetime reproductive success at the suggestive QTL on chromosome 26. This implies that sexual antagonism may play a role in maintaining genetic variation at this QTL. Our findings provide candidate regions for a classic avian life-history trait that will be useful for future studies examining the molecular and cellular function of, as well as evolutionary mechanisms operating at, these loci.